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Abstract 

The characterization of systems of differential equations admit- 
ting a superposition function allowing us to write the general solution 
in terms of any fundamental set of particular solutions is discussed. 
These systems are shown to be related with equations on a Lie group 
and with some connections in fiber bundles. We develop two methods 
for dealing with such systems: the generalized Wei-Norman method 
and the reduction method, which is very useful when particular so- 
lutions of the original problem are known. The theory is illustrated 
with some applications in both classical and quantum mechanics. 

1 Introduction 

Time evolution of many physical systems is described by non-autonomous 
systems of differential equations 

^ = X%x), ^ = l,...,n, (1) 

for instance, Hamilton equations or Lagrange equations when transformed 
to the first order case by doubling the number of degrees of freedom. 

The Theorem of existence and uniqueness of solution for such systems 
establishes that the initial conditions x(0) determine the future evolution. It 
is also well-known that for the simpler case of a homogeneous linear system 
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the general solution can be written as a linear combination of n independent 
particular solutions, . . . 

X = . . . , h,..., kn) = h H \-K X(n) , (2) 

and for each set of initial conditions, the coefficients can be determined. For 
an inhomogeneous linear system, the general solution can be written as an 
affine function of (n + 1) independent particular solutions: 

X = F{x(^i), ... , ki,... , k„) 

= + ki{x(^2) - + . . . + kn{x^n+i) - X(i)) . (3) 

Under a non-linear change of coordinates both system become non-linear 
ones. However, the fact that the general solution is expressible in terms of 
a set of particular solutions is maintained, but the superposition function is 
no longer linear or affine, respectively. 

The very existence of such examples of systems of differential equations 
admitting a superposition function suggests us an analysis of such systems. 
We are lead in this way to the problem of characterizing the systems of 
differential equations for which a superposition function, allowing to express 
the general solution in terms of m particular solutions, does exist. The 
solution of this problem is due to Lie [0]. Our aim here is to review the 
theory developed by Lie from a modern geometric viewpoint and to present 
different applications both in mathematics and physics. 

The paper is organized as follows: Section 2 present the main Theorem 
due to Lie and some simple examples are given in Section 3. In Section 4, 
after the introduction of some notation concerning ingredients of Lie group 
theory, the particular case in which the systems are defined in a Lie group 
G is analyzed, and we show how to relate them with a particular type of 
equations on a group. We also show that Lie systems in homogeneous spaces 
arise naturally associated with these systems in Lie groups. The theory is 
illustrated with a pair of examples which point out the universal character of 
the equation in the group. Section 5 is devoted to present a generalization to 
the general case of a method proposed by Wei and Norman for linear systems 
and the example of the affine group in one dimension is used to illustrate the 
theory. The relation of the problem at hand with the theory of connections 
is studied in Section 6: it is shown that Lie systems define horizontal curves 
with respect to a connection. The reduction method developed in section 7 
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corresponds to consider the action of the group of automorphisms of the prin- 
cipal bundle on the set of connections, transforming in this way the given 
problem in a simpler one. Some examples and references to different ap- 
plications of this reduction method are also given. The applications of the 
general theory to different problems in both Classical and Quantum Mechan- 
ics are indicated in Section 8 with an especial emphasis on time evolution of 
time-dependent Hamiltonian systems. The example of the time-dependent 
linear potential model has been explicitly developed in Section 9 in both the 
classical and the quantum case. 



2 Lie Theorem 

The characterization of non-autonomous systems (|I]) having the mentioned 
property that the general solution can be written as a function of m inde- 
pendent particular solutions and some constants determining each specific 
solution is due to Lie. The statement of the Theorem, which can be found 
in the book edited and revised by Scheffers [0, is as follows: 

Theorem 1 Given a non- autonomous system of n first order differential 
equations like Pj, a necessary and sufficient condition for the existence of a 
function F : M'^™"^^-' — » such that the general solution is 

X = F{X(^i), . . . ,X(^rny, ki,... , kn) , 

with I a = 1, . . . , m} being any set of particular solutions of the system 
and ki, . . . ,kn, n arbitrary constants, is that the system can be written as 

— = Z,{t)e\x) + ■■■ + Zr{t)C\x) , z = 1, . . . , n , (4) 

where Zi, . . . , Z^, are r functions depending only on t and a = 1, . . . ,r, 
are functions of x = {x^, . . . , x"), such that the r vector fields in R" given by 

^^n)_^ a = l,... ,r, (5) 

i=l 

close on a real finite- dimensional Lie algebra, i.e. the vector fields Y^""^ are 
linearly independent and there exist real numbers, /"^ y, such that 

r 

7=1 
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The number r satisfies r < mn. For a geometric proof, see 0]. 

From tlie geometric viewpoint, tlie system of first order differential equa- 
tions (|I]) provides tfie integral curves of the t-dependent vector field on an 
n-dimensional manifold M 

1=1 

in the same way as it happens for autonomous systems and true vector fields, 
and the t-dependent vector fields satisfying the hypothesis of the Theorem 
are those which can be written as a t-dependent linear combination of vector 
fields, 



x(x,t) = 5^z„(t)r(' 



with vector fields y-"-* closing on a finite-dimensional real Lie algebra. They 
will be called Lie (or even Lie-Scheffers) systems. Many of its applications in 
physics and mathematics have been developed by Winternitz and coworkers 

i-ii. 



3 Some examples 

We have mentioned in the Introduction two types of systems of differential 
equations whose general solution can be written as described by Theorem |l|: 
homogeneous linear systems like 

^ = X:^%Wx\ z = l,...,n, (7) 

i=i 

for which m = n and the (linear) superposition function is given by (0), and 
the inhomogeneous ones, 

— = ^A%(t)x^+5^(t) , ^ = l,...,n, (8) 

i=i 

for which m = n + 1 and the (affine) superposition function is (^). 



4 



In the first case, the hnear system can be considered as the one giving 
the integral curves of the t-dependent vector field 



X=Y.A^,{t)x^ — 



(9) 



which is a linear combination with time-dependent coefficients, 



of the vector fields 



dx^ 



Notice that 



kl\ 



dx^ ' dx^ 



i,j = 1,... ,n. 



5'^ x^ 



d 
dx^ 



5^^ x^ 



_d_ 

dx' 



(10) 



(11) 



I.e. 



[Xjj, Xki] — 5*' Xkj — 5^^ Xii 



(12) 



which means that the vector fields {Xij}, with i,j = 1,... ,n, appearing 
in the case of a homogeneous system, close on a n^- dimensional Lie algebra 
isomorphic to the Qi{n, M) algebra. It suffices to compare these commutation 
relations with those of the g[(n, M) algebra. The latter is generated by the 
matrices Eij with elements {Eij)i^i = SikSji, which satisfy 

[Eij, Em] = Sjk Eii — 6ii E^j . 

Therefore, in this homogeneous linear case r = n"^ and m = n, hence the 
inequality r < mn is actually an equality. 

For the case of the inhomogeneous system (^, the time-dependent vector 
field is 



n In 



(13) 
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which is a hnear combination with t-dependent coefficients, 

n n 
ij=l 1=1 

of the V? vector fields (|l5) and the n vector fields 

d 

X, = — , z = l,...,n. (15) 

Now, these last vector fields commute among themselves 
[Xi,Xk] = , Vi, A; = 1,... 

and 

[Xij^Xk] = -5kj Xi , \fi,j,k = l,...,n. 

Therefore, the Lie algebra generated by the vector fields {Xij,Xk \ i,j,k = 
1, . . . , ra} is isomorphic to the (n^ + n)-dimensional Lie algebra of the affine 
group. In this case r = + n and m = n + 1 and the equality r = mn also 
follows. 

Another remarkable example, with many applications in physics, is that 
of the Riccati equation, which corresponds to n = 1 [11], |12|: 

= C2{t) x\t) + Ci{t) x{t) + Co{t) . (16) 

In this case r = 3 and 

e{x) = i, e{x) = x, e{x)=x\ 

while 

Zi (t) =Co{t), ^2 (t ) = Ci (t) , Zs{t)=C2{t) . 

The equation (p!6| ) determines the integral curves of the t-dependent vector 
field 

x(x,t) = c2(t)F(=^) + ci(t)r(2) + co(t)r(^) , 

where the vector fields Y^^^ , Y^'^^ , and Y^^^ in the decomposition are given by 

F« = ^, Y(^^=x^, Y('^=x'^. (17) 
ox ox ox 
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It is quite easy to check that they close on the following three-dimensional 
real Lie algebra, 

[F(i),r(2)] = yd) , [r(i),r(3)] = 2Y^^^ , [y^^\y^'^] = r^^) , (ig) 

i.e. thes[(2,M) algebra. 

It can be shown that, for the Riccati equation, m = 3, and hence, as 
r = 3 the equality r = mn holds. The superposition function comes from 
the relation 

:^:^^ = fc, (19) 

X — X2 X^ — X2 

or, in other words |Tl|, 

Xi{Xs - X2) + k X2iXi - Xs) 

^ = — 7 ^ n ^ — • (20) 

(X3 - X2) + k [Xi - X3) 

In particular, the solutions Xi, X2 and X3 are obtained for k = 0, 00, and 1, 
respectively. 

We show next a last example of nonlinear superposition formula for a 
specific Lie system, in order to illustrate how complicated the explicit ex- 
pressions can become. For the sake of brevity, we give only the result. 

Consider the differential equation system 

dx{t) _ , , , /,x , , /,x / 2 „,2 



dt 

dyit) 



h^{t)+h2{t)x + h{t) (x'-y'), 
b2it)y + 2hit)xy, (21) 



dt 

which determines the integral curves of the t-dependent vector field 

X{x, t) = hit) F(^) + 62 (t) Y^^^ + bs{t) Y^'^ , 
where now Y^^\ Y^'^\ and Y^^^ are given by 

i^« = ^, y^'^=^^ + y^, Y('^ = {x'-y')^ + 2xy^. (22) 
ox ox oy ox oy 

These vector fields satisfy the commutation relations ([l8|) , hence the previous 
system is a Lie system with associated Lie algebra s[(2,M), and with r = 3, 
n = 2. The number of particular solutions needed is m = 3. In fact, suppose 
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we know a set of particular solutions = {xj, i = 1, 2, 3 of the system 
([21|). Then the general solution can be written as 

X = X(2), X(3), ki, k2) = 

N 

y = F2{X(^i), X(^2), 3^(3), h, k2) = ^ , (23) 

where 

iV^ = xi{(a;2 - x^y + (?/2 - Vsf} 

+ki{xl X3 + (x3 - X2)xl + {yi - y2fx^ 

-X2{xl + (yi - y^f) - xi((a;2 - 2:3)^ + (^2 - Z/s)^)} 

+ A;2{x3(?/2 - yi) + ^2 (1/1 - 1/3) + (1/3 - y2){x\ + (yi - t/2)(l/l " I/3))} 

+ + A;2)a;2{(a;i - X3)' + (1/1 - ^3)'} , 



= yi{{x2 - xsf + {y2 - y^f} 
+ki{xl{y'i - yi) - xliyi + ^2) + 2a;2(a;3yi - xiy^) 

+2x1X3^/2 - {xl + {yi + ?/2)(yi - 2/3))(z/2 - ^3)} 
+A;2{xi(x2 - X3) + X2 X3 + X3(y2 - yl) 

/2|2 2\| / 2 2,2 2m 
-X2(X3 + 1/3 - 1/1) + ^1(^3 - ^2 + Z/3 - ?/2)l 

+{kl + /c^)?/2{(xi - X3)' + {yi - yzf} , 



D = {x2- x^f + (1/2 - y^f 

-2/ci{(xi - X3)(x2 - X3) + {yi - yz){y2 - 2/3)} 
+2A;2{x3(y2 - yi) + X2{yi - y^) + xi{y^ - 2/2)} 
+{kl + A;2){(a;i - X3)' + {yi - ^3)'} , 

and fci, A;2 are two arbitrary real constants determining each particular so- 
lution. For example, the particular solutions {xi, yi}, {x2, 1/2} and {X3, y^} 
can be obtained by taking fci = ^2 = 0, the limit A;i ^ 00 (or k2 — * 00), and 
fci = 1, ^2 = 0, respectively. 

In particular, if we look for solutions of the system ( pi] ) with ?/ = 0, 
we recover, essentially, the Riccati equation ([T6|); likewise, the superposition 
formula (|23|) reduces to (^OD in such a particular case. 

For a more complete information about the explicit construction and use 



of superposition formulas, see, e.g., H]-[ll] and references therein. 
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4 Lie— SchefFers systems on Lie groups 



The most important example, which will be shown to give rise to many other 
related systems, occurs when M is a Lie group G and we consider vector fields 
Xa in G that are either left-invariant or right-invariant as corresponding 
either to the Lie algebra g of G or to the opposite algebra ||13|, |14|- 

Let us choose a basis {ai, . . . , ar} for the tangent space T^G at the neutral 
element e G G, and denote {di, . . . , 'dr} the corresponding dual basis oiT*G. 
In the following denotes the right-invariant vector field in G such that 
X^{e) = tta, i.e. 

and in an analogous way, X^ will denote the left-invariant vector field 

Similarly, 9^ and 9^ are the right- and left-invariant 1-forms in G determined 
by i.e. 

9^{g) = {Rg-.YM , 9^{g) = {Lg-^YM . 
If we consider the right-invariant Lie-Scheffers system on G 

r 

X{g,t) = -J2Ut)X^{g) , (24) 

a=l 

its integral curves will be determined by the system of differential equations 

r 

m = -Y.^.{t)X^i9{t)) . (25) 

= 1 

Applying {Rg{t)-'^)*g(t) to both sides of this equation we obtain the equivalent 
equation 

r 

{Rg{t)-^)*git){9{t)) = - ^ ba{t) a„ , (26) 

Q=l 

which we will write as well, with a slight abuse of notation, as 

r 

{gg-'m = -J2Ut)aa, (27) 

o=l 
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although reduces to { ^7\) only when G is a matrix group. This equation 
is right-invariant, and so, out of a solution g{t) of (^) with initial condition 
g{0) = e, the solution with initial conditions g{0) = go is given by g{t)go. 
This means that for the Lie-Scheffers system (^) on the Lie group G, m = 1. 

Of course, given a homomorphism of Lie groups F : G —>■ G', the right- 
invariant Lie-Scheffers system on G (^4]) produces a right-invariant Lie- 
Scheffers system on G', 

r 

X{g\t) = -Y,Ut) iF.X)2{g'), 

a=l 

where is the right-invariant vector field on G' which is F-related with 

the vector field X^. 

Let us consider a left action of a Lie group G, with Lie algebra g, on a 
manifold M, $ : G x M ^ M, and denote $g : M ^ G and $p : G G, 
where g E G, p & M, the maps defined by ^g{p) = ^p{g) = ^{g,p)- The 
fundamental vector field Xa associated to the element a of g is given by 

{XJ){p) = -^/($(exp(-to),p)) , / G G^{M) , 
at t=o 

where the minus sign has been introduced for X : g ^ X(M) to be a Lie 
algebra homomorphism, i.e. a M-linear map such that X[a_b] = [X^, Xf,]. Note 
that Xa{p) = $p^,e(— a). These vector fields are always complete. As an ex- 
ample, let us consider the left action action of G on itself by left translations, 
^{g,g') = g g' . The fundamental vector fields Xa are right invariant because 

{Xa){g) = $g*e(-a) = Rg.e{-a) = -{X^){g) , 

where X^ is the right-invariant vector field in G determined by its value at 
the neutral element (X^)(e) = a. 

Given two actions $i and $2 of a Lie group G on two differentiable 
manifolds Mi and M2, a map F : Mi M2 is said to be equivariant if 
F o $ig = o F. The remarkable property is that when G is connected, 
the map F : Mi — > M2 is equivariant if and only if for each a G TgG 
the corresponding fundamental vector fields in Mi and M2 are F-related 



Now, let if be a closed subgroup of G and consider the homogeneous 
space M = G/H. Then, G acts on M by X{g',gH) = {g'g)H. Moreover, G 
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can be seen as a principal bundle {G,t,G/H) over G/H, where r denotes 
the canonical projection. The important point is (see e.g. |]13|) that the map 
r : G — > G/H is equivariant, with respect to the left action of G on itself by 
left translations and the action X on G/H, and consequently, the fundamental 
vector fields corresponding to the two actions are r-related. Therefore, the 
right-invariant vector fields are r-projectable and the r-related vector 
fields in M are the fundamental vector fields —X^ = —Xa^ corresponding to 
the natural left action of G on M, T^gX^[g) = —Xa[gH). In this way we 
will have an associated Lie-Scheffers system on M: 

r 

X{x,t) =Y,ba{t)X^{x) , (28) 

a=l 

where x = gH, whose integral curve will be determined by 

r 

X = Y,ba{t)Xa{x) . (29) 



Thus, the solution of (^) starting from Xq will be x(t) = $((yf(t), Xq), with 
g{t) being the solution of ( pG] ) with g{0) = e. This is the main point: the 
knowledge of one particular solution of (p6|) allows us to obtain the general 
solution of (pHf). 

The converse property is true in the following sense: Given a Lie-Scheffers 
system in a manifold M defined by complete vector fields and with associated 
Lie algebra g, we can see these as fundamental vector fields relative to an 
action given by integrating the vector fields. Then, the restriction to an orbit 
will provide a homogeneous space of the above type. The choice of a point 
xo in the homogeneous space allows us to identify the homogeneous space M 
with G/H, where H is the stability group of xq. Different choices for xq will 
lead to conjugate subgroups [l^ . 



For instance, the vector fields appearing in ( piTj ) close on a Lie algebra but 
the third one is not complete on M. We can however consider the one-point 
compactification ofM, M = MU{oo}, and then the fiows of vector fields in 
([T7|) are, respectively. 



X 

X t-^ X + e , X ^ e X , x 



1 — xe 

and therefore they can be considered as the fundamental vector fields cor- 
responding to the action of S'L(2,M) on the completed real line M, given by 
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^/ . N ax + (3 .„ , 5 

^A,x) = — , lix^ — , 

7 a; + 7 

oo) = -, $ (A, -5/7) = oo, 

7 

when ^ = ( " ^ ) e ^L(2,M). (30) 

The stabihty group of 00 is the subgroup of matrices with 7 = 0, which is 
isomorphic to the affine group Ai in one dimensional space, while the stability 
group of is made up by the matrices with /? = 0, a group isomorphic to Ai- 
Indeed 

'5 0\ /O l\/a7\/0-l 



-1 a ) \-\ ) \ ^ b ) \ \ 

The remarkable fact is that equation p6| ) has a universal character. There 
will be many Lie-Scheffers systems associated with such an equation. It is 
enough to consider homogeneous spaces and the corresponding fundamental 
vector fields. In this way we will get a set of different systems corresponding 
to the same equation on the Lie group G. In particular, we can consider 
an action of G on a linear space given by a linear representation, and then 
the associated Lie systems are linear systems. Hence, we obtain a kind of 
linearization of the original problem 0. Therefore, the theory can be useful 
in the study of both classical and quantum problems. 

As an example we can consider both the Riccati equation 

and the linear system of first order differential equations 

d_( x\^( hit) boit) \fx 
dt\y ) V -Ht) -hit) J\y 

They are two different Lie systems for G = SL(2,M.) corresponding to the 
same equation 

99'' = -boit) Mo - 2 hit) Ml - hit) , 
with the matrices 



0/' ' 210-1/' ' V -1 



12 



being a basis of s[(2,M). They satisfy the commutation relations 



[Mo, Ml] = -Mo , [Mo, M2] = -2 Ml , [Mi, M^] = -M^ . 

As another illustrative example, we can consider the non-relativistic dy- 
namics of a spin 1/2 particle, when only the spinorial part is considered 



14| . The dynamics of such a particle in a time-dependent magnetic field is 



described by the so called Schrodinger-Pauli equation: 

^h'^ = H\iJ) = -^^B.S\iJ) , 

with n proportional to the Bohr magneton, B = {B^, B^, B^) the t-dependent 
magnetic field, and S'j = More explicitly. 



d_fipi\_[J^f iB^ iB^ + B^ \ f ipi 
dt[ ^2 J 2[ iB^ - B^ -iB^ j W: 



2 



(31) 



The matrices —ia^, —ia"^ and —ia^ generate the real Lie algebra of trace- 
less skew- Hermit ian 2x2 matrices, the Lie algebra of the group SU{2, C) 
and therefore of 5*0(3, M). 

As a consequence of the theory we have developed in this section, in 
order to find the general solution of the evolution equation (^Tj), it suffices 
to determine the curve R{t) in S'0(3,]R) starting from the identity map, 
R{0) = I and such that 



RR-^ = B^M3 + B^M2 + B^Mi 



where 





Ml = I -1 I , M2 = I 1, M3 




1 

Such a curve gives us the general solution for the dynamics 

\m) = Rmm , 

where R is an element in SU{2, C) corresponding to R. 
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5 The Wei— Norman method 



In order to solve directly the equation (pq ) we can use a method which 
is a generalization of the one proposed by Wei and Norman [jl5|, |l^] for 
finding the time evolution operator for a linear systems of type dU {t) /dt = 
H{t)U{t), with U{0) = I, see also However, as it will be mentioned in 
a later section, there exist alternative methods for solving (^) by reducing 
the problem to a simpler one. 

Both procedures are based on the following property ||13[: If g{t), gi(t) 
and g2{t) are differentiable curves in G such that g{t) = gi{t)g2{t), Wt G M, 
then, 

Rg{t)-^ *g(t){g{t)) = -Rgi(i)-i*gi(t)(fi'lW) + Ad(5(i(t)) {-RgaW"! *92(t) (5'2(i) ) } • 

(32) 

The generalization of this property to several factors is as follows. Let 
now g{t) be a curve in G which is given by the product of other / curves 
9{t) = 9i{t)g2{t) ■ ■ ■ gi{t) = Yl\=i9iit)- Then, denoting hs{t) = ULs+i 9t{t) , 
for s G {1, . . . , Z — 1}, and applying (|3^) to g(t) = gi(t) hi(t) we have 



Rg{t)-^ *g{t){9{t)) = Rgiit)-^ *g,(t){9l{t)) + M{gi{t)) ^Rh^(^t)-i ,hi{t){hl{t)) j ■ 

Simply iterating, and using that Ad{gg') = Ad{g)Ad{g') for all g, g' & G we 
obtain 

Rg{t)-^*g{t){9{t)) = i?g,(t)-l,g,(t)(^l(t)) + Ad(^l(t)) {Rg^it)-^ *g^{t){g2{t))} 

1=1 \j<i / 

,{t)-i *g,{t) 

where it has been taken go{t) = e for all t. 

The generalized Wei-Norman method consists on writing the solution 
g{t) of (p6D in terms of its second kind canonical coordinates w.r.t. a basis 
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{oi, . . . , ttr} of the Lie algebra g, for each value of t, i.e. 

r 

g{t) = JJexp(-t;„(t)aQ) = exp{-vi{t)ai) ■ ■ ■ exp{-Vr{t)c 



a=l 



and transforming the differential equation (pB]) into a differential equation 
system for the Va{t), with initial conditions Va{0) = for all a = 1, ... , r. 
The minus signs in the exponentials have been introduced for computational 
convenience. Then, we use the result (^3]), taking I = r = dimG and Qaif) = 
exY>{-Va{t)aa) for all a. Now, since Rg^{t)-^ *g^{t){ga(t)) = -Va{t)aa, we see 
that ( p3D reduces to 

Rg{t)-^ *git){9{t)) = -^Wa I JJ Ad(exp(-t;/3(t)a/3)) j a„ 

"=1 \/3<a / 

= - I n exp(-W/3(t)ad(a;3)) J , 

a=l \/3<a / 

where it has been used the identity Ad(exp(a)) = exp(ad(a)), for all a E Q. 
Substituting in equation (|26|) we obtain the fundamental expression of the 
Wei-Norman method 

I J]^ exp(-t;/3(t)ad(a/3)) J = ^&„(t)a„ , (34) 

a=l \l3<a / a=l 

with Va{0) = 0, a = 1, . . . , r. The resulting differential equation system for 
the functions Va(t) is integrable by quadratures if the Lie algebra is solvable 



15| , |T6[, and in particular, for nilpotent Lie algebras. 

As a simple but illustrative example we can consider the affine group in 
one dimension, Ai, i.e. the set of transformations of the real line 

X = ai X + ao , (35) 

with ai 7^ and ao being real numbers. The group composition law is 

(«[), a[) * (ao, tti) = (tto + '^1 '^0, a'l «i) • 
Denoting by {xo,Xi) the coordinate system in Ai given by 
a;o(ao,ai) = ao, Xi(ao, ai) = ai , 
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we see that a basis of right-invariant vector fields in Ai is given by 



dxo ' dxo dxi 

while the corresponding basis of left-invariant vector fields is given by 

therefore the defining relations are 
Then, 

ad(ao)ao = , ad(ao)ai = — , 
and ii g = exp(— uoao) exp(— uiai), equation ( ^4t) becomes in this case 

Mo «o + + % ^o) =^0^0 + ^1 Oi 5 
so we obtain the system 

Mo = &o - h Uq, Ml = 6i , (36) 

with the initial conditions uo{0) = ui{0) = 0. Note that the first equation is 
nothing but an inhomogeneous linear equation. The explicit solution can be 
obtained through two quadratures: 

(t) = e-Iodt'Mt') f dt'bo{t')e^odt"hit") ^ ^ fdt'hit') , 



In a similar way, if we consider instead g = exp(— wiOi) exp(— foOo), and 
we take into account that 

ad(ai)ao = clq ? ad(ai)ai = , 

then we will find 

vi ai + vq exp(— t>iad(ai))ao = Viai + vq e"*'^ oq = 60^0 + oi , 
yielding the system 

i'o = bo e''! , vi = bi , (37) 



16 



also with the initial conditions fo(0) = fi(0) = 0. The system (|T7D, with 
such initial conditions, can be easily integrated by two quadratures: 



When we consider the previous action (^) on the real line we get as 
fundamental vector fields Xq = — ^ and Xi = — thus the Lie system in 
M which corresponds to g g^^ = —bo oq — 6i ai is x = —bo — bi x. Our theory 
gives us the formula for the explicit general solution of such inhomogeneous 
linear differential equation, by making use of 

x{t) = <l>(5((t),Xo) = $(exp(-Mo(t)ao),$(exp(-Mi(t)ai),Xo)) 
= e-"^Wa;o-Mo(t) , 

namely 

x{t) = e-/o jxo - j\t'bo{t') e/o'^*"^i(*")| . (38) 

Likewise, the same solution can be obtained from the second factorization: 

a;(t) =$(exp(-i;i(t)ai),$(exp(-i;o(t)ao),a;o)) = e-"^(*)(2;o-t^oW) , (39) 
which clearly gives the same result. 



6 Connections and Lie systems 

If G is a connected Lie group, the set of curves 

7 : M^G 
t ^ g{t) 

is also a group when the following composition law is considered: 

72 *7i : t g2{t)gi{t) . 

Given a curve 7 in G, g{t), such that g{^) = e, then g{t) = g(t) go is 
another curve in G starting from go (it is said to be right translated of 7 by 
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Qq) and similarly g{t) = gog{t) is also a curve in G starting from qq (called 
left-translated from 7 by go). 

Now, given the curve 7, we have a vector field along 7 given by the tangent 
vector g{t), and then, translating these tangent vectors to the neutral element 



by Rg-^(t)*g{t) we obtain a curve in T^G like in (|27|). The curve g(t) = g{t) go, 
right-translated by go of 7, gives rise to the same equation. 

But we can consider such equation ( ^7]) as an equation for the curve g{t) in 
G determined by the curve a(t) = — Yla=i ^aif) in T^G. This equation is 
right-invariant in the sense that if g{t) is a solution such that (?(0) = e, then, 
for each go G G, g{t) = g{t) go is a new solution, now such that g{Q) = go. 

We remark that if G is a Lie group, tt2 '. P = G x ^ defines a 
principal G-bundle. The right action of G on P is given by 

<iliig',t),g) = ^il,iig',t)) = ig'g,t) , 

i.e. '^g = RgX idiR. 

Giving a connection in P is equivalent to give a curve in G, for instance, 
one such that g{0) = e. It is also well-known that each global section provides 
a different trivialization of the principal bundle P. The given curve furnishes 
a section for 712, cr(t) = {g{t),t), and a family of sections right-translated 
from such a section, 

{a\t) = ^{a{t),go)\goeG} . 

The tangent vectors to such family of sections span the horizontal spaces 
in each point. More specifically, horizontal and vertical spaces in a point of 
P are given by: 

yP(,„,) = ((Xf((?o),0)) , 
HP^g,,) = {iRgUm9~\t))A)) ■ 

Note that 

The choice of the connection given by 7 amounts to choose a basis of the 
tangent space at the point {go, t) as follows, {X^{go), d/dt+Rg^^eigit) 9^^(t))}: 
while the dual basis is made up by {Oa{9o)—Ta(t) dt, dt] where the coefficients 
Ta{t) are determined by 

Rg,.e{g{t)9'\t)) =Y,rp{t)X^{go) , 

/3 
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i.e. Ta(t) = {0^{go), Rgg^e{g(t) 9 ^(t)))- Therefore, the vertical projector 
associated to the connection is 

^i9o,t) = J2^?(9o) ® (Oj^igo) - rpit) dt) = idT,,G - {RgoUm 9-\t))) dt . 



It is also well-known that when a left action ^ : G x M ^ M of G on M 
is considered, there exists an associated bundle E with base M and typical 
fiber M. The total space of such bundle is the set of orbits of the right action 
of G on P X M, 

{u,x)g = {'^{u,g),<l>{g-\x)) , 

being the projection 7!-e[u,x] = Tr2{u), where [u,x] denotes the equivalence 
class of {u,x) E P X M and u is of the form {g',t). A connection in the 
principal bundle translates into a connection in the associated bundle E, 
and so the horizontal curves will then be [7(t), x], where 7(t) is an horizontal 
curve in P. More exphcitly, as the curves 7 are of the form 7(t) = {g{t)go, t), 
we find that the horizontal curves in the associated bundle are 

[{g{t)go,t),xo] = [^{{e,t),g{t)go),Xo] = [{e,t),<!>{g{t)go,Xo)] , 

and consequently, 

[{g{t)go,t),xo] = [(e,t),<l>(^(t),$(^o,a;o))] • 

Since the principal bundle is trivial, E is equivalent to a product. When 
$ is transitive, E = M x M., where [(e,t),a;] corresponds to {x,t) and with 
this identification, the horizontal curve here considered corresponds to the 
integral curve starting from the point ^{go,xo) of the associated Lie system 
in M with respect to the action of G on M given by $. 

Of course the simplest case is when a linear representation of G on the 
vector space V is considered, the associated bundle being then a vector bundle 
and the corresponding Lie system being a linear system. This means that a 
system as in (|^) can be seen as defining the horizontal curves corresponding 
to a connection in an associated vector bundle. The fact that linear systems, 
as Schrodinger equations, could be thought of as defining horizontal curves of 



a connection were considered several years ago [17| and it has been suggested 



recently by looking at the transformation properties of the equation under 



certain gauge changes [IS 
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7 The reduction method 



Given an equation on a Lie group 

r 

g{t)9{ty' = <t) = -J2ba{t)aa G T,G , (40) 



a=l 



with g{0) = e ^ G, it may happen that the only non- vanishing coefficients 
are those corresponding to a subalgebra f) of g. Then the equation reduces 
to a simpler equation on a subgroup, involving less coordinate functions in 
the Wei-Norman method. 

On the other hand, we know that such an equation pO|) can be seen as a 
connection in a principal bundle, and it is also well-known that the group of 
automorphisms of the principal bundle acts on the set of connections. The 
automorphims of the bundle we are considering are given by curves g'{t) in 
the group G, and so the group of curves in G defines an action on the set 
of connections and therefore on the Lie systems on the group. We can take 
advantage of such an action for transforming a given Lie system in another 
simpler one. 

Now, let us choose a curve g'(t) in the group G, corresponding to a given 
automorphism, and define the curve g(t) by g(t) = g'(t)g(t), where g(t) is 
the previous solution of (P^). The new curve in G, g{t), determines a new 
connection and therefore a new Lie system. 

Indeed, from 



Rg{t)-^*9it)(9{t)) = Rg'-Ht)*g'(i){g{t)) - ^ ha{t)Kd{g'{t))ao, , (41) 

which is an equation similar to (^) but with different right hand side. There- 
fore, the aim is to choose the curve g'{t) appropriately, i.e. in such a way 
that the new equation be simpler. For instance, we can choose a subgroup 



H and look for a choice of g'{t) such that the right hand side of (|T) lies in 
T(,H, and hence 'g{t) G H for all t. 

Now, suppose we consider a transitive action $ : G x M — ^ M of G 
on a homogeneous space M, which can be identified with the set G/H of 
left-cosets, by choosing a fixed point Xq: H is then the stability subgroup 
of xq. The horizontal curves starting from the point xq associated to both 
connections are related by 

x{t) = $(^(t),xo) = <l>{g'{t)g{t),xo) = ^g'{t),<l>{g{t),xo)) = ^g'{t),x{t)) . 
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Therefore, the action of the group of curves in G on the set of connections 
translates to the homogeneous space and gives an action on the corresponding 
set of associated Lie systems. More exphcitly, if we consider the automor- 
phism defined by g'{t), the Lie system (^) transforms into a new one (see 

H) 

r 

^ = Y,ba{t)M^) , (42) 

a=l 

in which 

b = Adig'it))b{t)+g'g'-' . 

The important result proved in [T^ is that the knowledge of a particular 
solution of the associated Lie system to ( p6D in G/H allows us to reduce 
the problem to one in the subgroup H. For any choice of the curve g'{t) we 
can consider a curve x'{t) defined in the homogeneous space G/H by g'{t) as 
follows: x'{t) = T{g'-^{t)g{t)) = g'~^{t)H. Then, if g'{t) is chosen such that 
the curve x'{t) is a solution of the associated system, then the automorphism 
defined by g'{t) transforms the original problem into one in the subgroup H: 

Theorem 2 Every integral curve of the time- dependent vector field on the 
group G, given by the right hand side of can be written in the form 

g{t) = gi{t)h{t), where gi{t) is a curve projecting onto a solution Xi{t) of 
an equation of type (j^^ for the natural left action of G on the homogeneous 
space G/H , and h(t) is a solution of a type ^dj ) equation but for the subgroup 
H , given explicitly by 

{hh-'m = -Ad{g^\t)) {^K{t)a^ + {gigi'){t)^ e T,H . 

As an example, we can consider once again the affine group in one di- 
mension, Ai, of the preceding section. We can choose first the Lie subgroup 
Hq = {(ao,l) I ao G M} and consider the corresponding one-dimensional 
homogeneous space Ai/Hq. Its points can be characterized hj y = Xi, with 
Xi 7^ 0. In this coordinate system for Ai/Hq the fundamental vector fields are 
Xo = and Xi = y^. The Lie system associated to g g~^ = —b^ — bi ai is 
y = biy and according to the result of the preceding Theorem, once we know 
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a solution of this last homogeneous linear equation, we can carry out the 
reduction procedure. More explicitly, when we know a solution oi y = biy, 
the change of variable x = y( will simplify the equation x = bo + bi x to one 
on the subgroup Hq, ( = b^y'^. 

If we consider instead the Lie subgroup Hi = {(0,ai) | ai G M — {0}}, 
then the elements of the one-dimensional homogeneous space Ai/Hi can be 
characterized hj z = Xq. The expression of the fundamental vector fields in 
this coordinate system are Xq = and Xi = z-^. Then, as soon as we know 
a solution of i = 612; + 60, namely a particular solution of the inhomogeneous 
equation, we can reduce the problem of finding the general solution of i; = 
bo+bi X to solving an equation on Hi, which is a homogeneous linear equation. 
This procedure corresponds to the change of variables x = z + (, which leads 
to the reduced equation ( = bi(. 

Therefore the two methods usually found in textbooks for solving the 
inhomogeneous linear differential equation appear here as particular cases of 
a more general methodology for reduction of differential equation systems to 
simpler ones. 

In this way, the last method can be generalized when one considers an 
inhomogeneous linear system, whose associated group is the corresponding 
affine group. Given a particular solution, the problem is reduced to another 
one on its stabilizer, i.e. the group GL{n,M.), or, in other words, to a homo- 
geneous linear system. 

As another example, we consider the Riccati equation, which has been 
shown to be an example of Lie system, corresponding to the left action (|30D 
of SL{2,M.) on the (compactified) real line M by homographies, see e.g., 

|T2|. The action of the group of automorphisms of the principal bundle 
translates to the space R into an action of the group of curves in SL{2, R) 
on the set of Riccati equations. This action was used in [1^ for studying the 
integrability properties of the Riccati equation. The stabilizer of the point 
at the infinity is the affine group in one dimension Ai. Therefore, if we know 
a particular solution, xi, of the Riccati equation, the problem reduces to one 
on ^1, i.e. a inhomogeneous linear equation, by means of the well-known 
change of variable x = xi + z. Had we chosen the origin G R as the initial 
point, the stabilizer (isomorphic to ^1) would be generated by dilations and 
cotranslations. This corresponds to a new reduction of the Riccati equation 
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by means of the change of variable 

, X XXi 



X 



-— + 1 Xi- X 



which transforms x = Co(t) + Ci{t) x + C2(t) into 

x=-=^-^ + c,it)jx+coit), 
with associated group ^i. For more details, see [|12| and [|13] 



Now, suppose we know not only one but two different particular solu- 
tions, a;(i) and X(^2), of a Lie system in a homogeneous space. They will be 
determined by the choice of initial conditions which provide different presen- 
tations of the homogeneous space as G/Hi and G/H2, respectively, where 
Hi is the stability subgroup of a;(j)(0). Using the result of Theorem 2, with 
gi(t) being a lifting to G of both curves and a;(2)(t), we will get an 

equation like the one in the Theorem but where the right hand side will be 
in the intersection T^Hi fl T^H2, and therefore the Lie system is reduced to 
one on the subgroup Hi n The example of the Riccati equation was 
explicitly considered in |jl2[, where it was also shown that the knowledge of a 



third solution reduces the problem to a trivial equation a; = 0, and therefore 
giving rise in this way to the superposition function (pO|). 



8 Some applications in Classical and Quan- 
tum Mechanics 

Non-autonomous linear systems and Riccati equations are examples of Lie 
systems that appear very often in Physics. For instance, linear systems 
appear in the time evolution of time-dependent harmonic oscillators and 
the latter is a condition for the super-potential W in the factorization of a 
typical quantum Hamiltonian H = —(P/dx'^ + V{x) a.s H — e = {—d/dx + 



W){d/dx + W), where e is a constant (see, e.g., [|T^, 0), and it plays a 



relevant role in the search for the so-called Shape Invariant potentials (see 



1^ , pot P^ , |23| , p4|). As we have pointed out in preceding sections, the 
Riccati equation may appear each time that the group SL{2,'R) plays a 
role, and because of the isomorphism of the Lie algebras of SL{2,'R) and 
the linear symplectic group in two dimensions, it will be useful in the linear 
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approximation of symplectic transformations and the theory of aberrations 
in optics [pSD - 

However, the Riccati equation is particularly important because it ap- 
pears as a consequence of Lie reduction theory when taking into account 
that dilations are symmetries of linear second order differential equations 



TT|. Actually, the homogeneous linear second-order differential equation 

|£ + ,M| + eW. = 0, (43) 

admits as an infinitesimal symmetry the vector field X = z d / dz generating 
dilations in the variable z, which is defined for z 7^ 0. According to Lie 
theory we should change the coordinate z to a new one, u = (p{z), such that 
X = d/du. This change is determined by the equation Xu = 1, which leads 
to u = log \z\, i.e. 1^1 = e". In both cases of regions with 2; > or z < we 
have 

dz du ^ d'^z / du\ ^ d'^u 

dx dx ' dx"^ \dx ) dx"^ ' 



so the equation ( ^3]) becomes 



d^u , , .du ( du\ ^ , , 
+ ^\-r] +c(a;) = 



dx^ dx \dx ) 

and the order can be lowered by introducing the new variable w = du/dx. 
We arrive to the following Riccati equation for w 

^ = —w"^ — h{x)w — c{x) . (44) 
dx 

Notice that w = z~'^dz/dx, and that this relation together with (0) is 
equivalent to the original second order equation. In the particular case of 
the one-dimensional time-independent Schrodinger equation 

'^''^ + (V(x)-e)0 = O, 



dx^ 

the reduced Riccati equation for W = (j)~^d(j)/dx is 

W' = -W^ + {V{x) - e) , (45) 
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which is the equation that W must satisfy in the previously mentioned fac- 
torization of if = —(P/dx^ + V{x). 

Equations of type P5| ) are particular cases of Riccati equations. We have 
shown in a preceding section that it is possible to act with the group of curves 
in SL[2, M) on the set of Riccati equations. That means that a given Riccati 
equation can be transformed in other related equations (for more explicit 
details see, e.g., [Q). Therefore, when using curves in 5'L(2,R) preserving 
the form of a given equation like (0), we are transforming the spectral 
problem for a given Hamiltonian into that of another one. This method is 
carefully explained in |2^, where explicit examples of the usefulness of the 



theory are given. 

In a typical problem of Classical Mechanics we are dealing with Hamil- 
tonian vector fields in a symplectic manifold (M, Q) , and then we should 
consider the case in which the vector fields arising in the expression of the 
t-dependent vector field describing a Lie system, are Hamiltonian vector 
fields closing on a real finite-dimensional Lie algebra. These vector fields 
correspond to a symplectic action of the group G on the symplectic man- 
ifold (M, Q) . The Hamiltonian functions of such vector fields, defined by 
i{Xa)^ = —dha, however, do not close on the same Lie algebra when the 
Poisson bracket is considered, but we can only say that 



d {{ha,h,3} - h[x^,Xp]) = 



and therefore they span a Lie algebra extension of the original one. 

The situation in Quantum Mechanics is quite similar. It is well-known 
that the separable complex Hilbert space of states 7i can be seen as a real 



manifold admitting a global chart [^. The Abelian translation group allows 



us to identify the tangent space T^H at any point (f) G H with H itself, the 
isomorphism being obtained by associating with ip eH the vector ip G T^H 
given by: 

^/(0):= + ^ ' 

for any / G C°°(7^). 

The symplectic 2-form Q is given by 

with ("I-) denoting the Hilbert inner product on Ti. 
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Through the identification of 7i with T^Ti a continuous vector field is just 
a continuous map A: 7i ^ Ti; therefore a hnear operator A on 7i is a special 
kind of vector field. 

Given a smooth function a: 7i — > M, its differential da^ at G 7i is an 
element of the (real) dual H' given by: 



Now, as it was pointed out in [27| the skew-Hermitian linear operators 



in Ti, define Hamiltonian vector fields, the Hamiltonian function of —i A for 
a self-adjoint operator A being a(0) = | (0,^40). The Schrodinger equation 
plays the role of Hamilton equations because it determines the integral curves 
of the vector field —i H . 

Now, Lie system theory applies to the case in which a t-dependent Hamil- 
tonian can be written as a linear combination with t-dependent coefficients 
of Hamiltonians Hi closing on, under the commutator bracket, a real finite- 
dimensional Lie algebra. The remarkable point, however, is that this Lie 
algebra does not necessarily coincide with the corresponding classical one, 
but it is a Lie algebra extension. An example will be given in next section. 



9 An example: classical and quantum time- 
dependent linear potential 

The linear potential model, with many applications in physics, has recently 
been studied by Guedes We can use this problem in order to illustrate 
the possible applications of the theory by means of a simple example. Let us 
consider the classical system described by a classical Hamiltonian 

and the corresponding quantum Hamiltonian 

describing, for instance when f{t) = q Eq+q E cos ut, the motion of a particle 
of electric charge q and mass m driven by a monochromatic electric field. Eq 
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is the strength of the constant confining electric field and E that of the 
time-dependent electric field that drives the system with a frequency ujII-k. 
Instead of using the Lewis and Riesenfeld invariant method as it was 
done in p8|, we will study simultaneously the classical and the quantum 



problem by reduction of both problems to similar equations and using the 
Wei-Norman method to solve such an equation. As it will be shown, the 
only difference is that the Lie algebra arising in the quantum problem is not 
the same as in the classical one, but a central extension. 
The classical Hamilton equations of motion are 

V 
m 

P = -fit), (46) 

and therefore, the motion is given by 

, 1 /"t ft' 

x{t) = xo + — -- dt' fit") dt" , (47) 
•m m Jq Jq 

Pit) = po- [ fit')dt' (48) 
The t-dependent vector field describing the time evolution is 

m ox op 
This vector field can be written as a linear combination 

X = -X^-fit)X,, 
m 

with 

Xi=p — , X2 = — , 
ox Op 

being two vector fields closing a three dimensional Lie algebra with 

ox 

isomorphic to the Heisenberg algebra, namely, 

[Xi, X2] = -X3 , [Xi, X3] = , [X2, X3] = . (49) 
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The flow of these vector flelds is given, respectively, by 



01 (t, {xo,po)) = {xo+pot,po) , 
02(i, ixo,Po)) = ixo,Po + t) , 
03 (i:, (xo,Po)) = ixo + t,po) . 

In other words, this corresponds to the action of the Lie group of upper 
triangular 3x3 matrices on , 





It is to be remarked that the three vector flelds Xi, X2 and X3 are 
Hamiltonian vector flelds with respect to the usual symplectic structure, 
fl = dxAdp, the corresponding Hamiltonian functions hi such that i{Xi)Q = 
—dhi being 

hi = -—, h2 = x, h3 = -p, 

therefore 

{h,,h2} = -h, {h,,h^} = 0, {/l2,/l3} = -l, (50) 

which close on a four- dimensional Lie algebra with /i4 = 1, that is, a central 
extension of that given by (^Q]). Let {oi, 02, 03} be a basis of the Lie algebra 
with non- vanishing deflning relations [oi, 02] = —03. Then, the corresponding 
equation in the group ( P7| ) becomes in this case 

gg~^ = -— ai + f{t) . 
m 

Now, choosing the factorization g = exp(— ^303) exp(— ^202) exp(— uiai) 
and using the Wei-Norman formula (|3^) we will arrive to the system of 
differential equations 

m 

together with the initial conditions 

u,{0) = U2{0) = us{0) = , 
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with solution 







)i 









Ml = — , U2 = - f fit') dt' , U3 = -— f dt' f fit") dt" . 
m Jo mJo Jo 

Therefore the motion will be given by 

1 ^ -iiodt'i^'fxndt" 

1 -Iofit')dt' 

1 

which reproduces (^8]). So, obviously we will recover the constant of motion 
given in p8|, Ii = p{t) + /(f) dt', together with the other one I2 = x{t) — 

i {pit) + /o fit') dt') t + ^ dt' fit") dt". 

As far as the quantum problem is concerned, also studied in a very recent 
paper [PU|, notice that the quantum Hamiltonian Hg may be written as a 
sum 

Hg = -H,-fit)H2, 
m 

with 

p2 

H\ = — , Ho = — X , 
2 

and —i Hi and —i H2 close on a four- dimensional Lie algebra with —i H3 = 
—iP, and —iH^ = il, isomorphic to that of (^), which is an extension of 
the Heisenberg Lie algebra (^9]), 



[-1 Hi, -I H2] = -iH^, [-1 Hi, -I H,]=0, [-1 H2, -I H,] = -tH,. 

The Schrodinger equation given by the Hamiltonian Hq is like that of a Lie 
system. Note that this Hamiltonian is time-dependent and that such systems 
are seldom studied, because it is generally difficult to find the time evolution 
of such systems. However, this system is a Lie system and therefore we can 
find the time-evolution operator by applying the reduction of the problem to 
an equation on the Lie group and using the Wei-Norman method. 

Let {ai, a2, 03, 04} be a basis of the Lie algebra with non-vanishing defin- 
ing relations [01,02] = 03 and [02,03] = 04. The equation (^7|) in the group 
to be considered is now 

gg~^ = -— oi + fit) 02 . 
m 
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Using the factorization g = exp(— M4 04) exp{—us 03) exp(— 02) exp{—Ui Oi) 
the Wei-Norman method provides the following equations: 

«l = ^, U2 = -f{t), 

1 2 

Us + U2Ul = 0, U4 + U3U2 - -U2U1 ^ , 

and written in normal form 

^1 = ^ , ^2^ -f{t) , 

U3^-^U2, f(t)u3+ —ul, 

together with the initial conditions Ui{0) — ^2(0) — ^3(0) = ^4(0) — 0, 
whose solution is 

Ul{t) = - , U2{t) = - f fit') dt' , Usit) = - f dt' f fit") dt" , 

m Jo rnJo Jo 



and 



U4 

m 

These functions provide the explicit form of the time-evolution operator: 

C/(i,0) = exp(-W4(i))exp(ra3(t)P)exp(-ra2(t)^)exp(wi(t)P72) ■ 

Notwithstanding, in order to find the expression of the wave-function in 
a simple way, it is advantageous to use the factorization 

g — exp(— 1;4 04) exp(— 1;2 02) exp(— f 3 03) exp(— Ui Oi) . 

In such a case, the Wei-Norman method gives the system 

^'1 = ^ , ^2 = -fit) , 
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jointly with the initial conditions fi(0) = f2(0) = ^3(0) = f4(0) = 0. The 
solution is 



Then, applying the evolution operator onto the initial wave-function ipip, 0), 
which is assumed to be written in momentum representation, we have 



ij{p,t) = f/(t,0)V'(p,0) 

= expi-iv^it)) exp{-iv2it)X) exp{iv3{t)P) exp{ivi{t)P'^ /2)^/j{p, 0) 

= exp{-ivi{t)) exp{-iv2it)X)e'^'"^'^P+^'^'^P'/^'^iP{p, 0) 

= exp(-it;4(t))e*^"^^^*)^^+"'^*^^+"'^*^^^+"'^*^^'/'V(p + ^2(t), 0) , 



where the functions Vi(t) are given by ([^), (|^ ) and (p^D, respectively. 
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